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Recent works have explored the properties of Levy flights with resetting in one-dimensional do¬ 
mains and have reported the existence of phase transitions in the phase space of parameters which 
minimizes the Mean First Passage Time (MFPT) through the origin [Phys. Rev. Lett. 113, 220602 
(2014)]. Here we show how actually an interesting dynamics, including also phase transitions for the 
minimization of the MFPT, can also be obtained without invoking the use of Levy statistics but for 
the simpler case of random walks with exponentially distributed flights of constant speed. We ex¬ 
plore this dynamics both in the case of finite and infinite domains, and for different implementations 
of the resetting mechanism to show that different ways to introduce resetting consistently lead to a 
quite similar dynamics. The use of exponential flights has the strong advantage that exact solutions 
can be obtained easily for the MFPT through the origin, so a complete analytical characterization of 
the system dynamics can be provided. Furthermore, we discuss in detail how the phase transitions 
observed in random walks with resetting are closely related to several ideas recurrently used in the 
field of random search theory, in particular to other mechanisms proposed to understand random 
search in space as mortal random-walks or multi-scale random-walks. As a whole we corroborate 
that one of the essential ingredients behind MFPT minimization lies in the combination of multiple 
movement scales (whatever its origin). 


Resetting can be defined as the interruption of a process in such a way that it is brought back instantaneously from 
its present state to another fixed state (presumably the initial one) and allowed to evolve once again from there. In 
the last years there has been a growing interest for studying classical and well-known stochastic tran^ort processes 


ran^o 

(e.g. Brownian motion, Continuous-Time Random Walks, etc) when subject to stochastic resetting [ij-y ■ This idea of 
allowing the particles to start anew its process after some random time has been claimed to have potential applicability 
to the study of many different systems. For instance, animal foraging is often constrained by the existence of a homing 
dynamics that makes individuals return from time to time to its nest in , leading to the ecological concept of central 
place foraging. Alternatively, the tendency to revisit well known places in animals can be more complex and be 
driven by sophisticated memory mechanisms [^, a fact which is of particular interest since similar ideas have been 
used to explain human mobility patterns, too [lOl 111 ]. These behavioral mechanisms, driven by higher cognitive 
processes, require an extension of the concept of resetting whose properties have just began being understood [l2Lll3l|. 
In a different context, the transport of information packets through both wired or wireless networks is often subject 
to transmission losses (errors, overflow, ...), which may require the transmission of new packets after some time (a 
mechanism which also coincides with resetting under proper assumptions) for increasing the probability of reaching 
a target whose location is unknown h4|. Similar mechanisms could apply too to beam propagation through optical 
fiber with random inhomogeneities [3| and many other transport processes on graphs and complex media. 

Beyond this practical interest, processes with resetting are attractive from a theoretical perspective as they can 
induce dramatic changes in the system dynamics if compared to their resetting-absent counterparts. So, resetting 
facilitates the emergence of nonequilibrium stationary states in different scenarios as in free diffusion 0, [T^ , in motion 
subject to potential landscapes 


can induce heavy-tailed decays 


.equilibrium stationary states in qitterent scenarios as in 

__ [l^, or in coagulation-diffusion process es IitI . Also, it_^ _ 

m the probability distribution of models of transport and/or growth [ll.ll8l|. On more mathematical grounds, resetting 
processes have been found to have interesting connections to some variants of the multi-armed bandit problem [l9l.[^. 
which may contribute to extend their potential range of interest. 

The focus of the present work is in the use of stochastic resetting as a possible mechanism for optimizing random 
spatial searches. The relevance of this problem has been already recognized and discussed in some of the aforemen¬ 
tioned works (see, e.g., 0,ii) but a more specific discussion, addressing the connections of resetting to well known 
concepts and results from the literature on random search theory, is still lacking. We aim thus at covering this existing 
gap. 

First of all, let us note that search efficiency can be measured in many different ways attending to time, energetic, 
or other biological or physical constraints [2lj. In agreement with most statistical physics approaches to the field we 
identify here efficiency with the time required to reach the target, so we will focus on the Mean First Passage Time 
(MFPT) of the random walker through the target location as the main magnitude of interest. Previous works have 
already determined the existence of an optimal value of the resetting rate r which minimizes the MFPT of Brownian 
walkers through the origin in a semi-infinite media 0,i- The extension of this idea to the case of particles moving 
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according to Levy Flights has revealed the existence of a hrst-order transition in the optimal values of the phase space 
(r, p.) (with fj, is the Levy index) for which the MFPT becomes a minimum. So, for small values of the initial position 
Xq intermediate optimal values = (0.25,1.2893...) are found. This r^ime extends up to a critical value (xg)^,. 

such that for xg > (a^o)cr optimal combination becomes (0.22145..., 0) [J]. Here we will revisit this problem and 
will show how actually hrst-order phase transitions (although of a different nature) do also appear also for simpler 
random-walks, in particular for walks with exponentially distributed flights in hnite domains. Furthermore, we will 
try to establish a connection between these results and other models as mortal random-walks and multi-scale random 
walks which have been used previously in the random search literature. As a whole, we will show the existence of 
a common force behind all these motion mechanisms driving MFPT minimization. As we argue here, this force is 
strongly related to the compromise between efficiently covering nearby regions (exploitation) and the exploration of 
new areas. From this perspective, resetting must be viewed as a mechanism that promotes exploitation as it allows 
the possibility to revisit regions which may have been initially missed. This is exactly the same idea behind Levy or 
multi-scale search strategies [nil. So that, in the forthcoming Sections we will try to understand the reach of this 
analogy as a way to obtain general principles driving optimization of search processes. 


FIRST-PASSAGE PROPERTIES OF RANDOM-WALKS WITH RESETTING 


We will consider for the moment random walkers moving in a one-dimensional infinite domain (which under proper 
assumptions provide a convenient representation of many realistic search processes [illii) , and will focus on the 


role that stochastic resetting plays on the MFPT through the origin, given the walker starts from an arbitrary initial 
position xq. Note in advance that resetting is not expected to be a convenient mechanism for searching unless the 
target is relatively close to ccq; otherwise resetting would just prevent the walker from reaching new unexplored areas. 
So that, the use of an infinite domain with resetting can effectively serve to describe the search for nearby targets in 
finite domains provided that the domain size is much larger than the typical distance covered before resetting. 

In contrast with the cases studied in previous works, we will rather consider here that walkers move according to an 
isotropic velocity model. So, the flights will not be done instantaneously after some waiting time but the transition 
occurs progressively at a fixed velocity -bwo or —vg, each with probability 1/2. We do so because we consider that a 
velocity model is probably a more realistic choice for most search applications; anyway, we stress that the extension 
of our results to the case of instantaneous jumps separated by waiting times would be straightforward. 

We will compare the results for two different resetting mechanisms for the sake of completeness. The first one consists 
of giving the walker the possibility to reset its position to xg whenever a single flight is completed. Accordingly, we 
will denote by Xi, A 2 ,... the successive positions of the particle after the first, second,... event (where each event can 
be either a flight or a reset). Thus the position of the particle after the (i + l)-th event is chosen according to the rule 


f Xg, with prob. r if the i-th event was not a reset 
\Xi±vgUi, otherwise 


Here the flight durations C/i’s are independent and positive random times determined by the probability distribution 
function (pdf) <pit), and the symbol ± denotes that the increment in Xi can be positive or negative with the same 
probability. In the following, we will term this resetting mechanism as being subordinated to flights since the statistics 
of the flights determines in part the rate at which resetting will occur. Note that this mechanism is slightly different 
to that proposed in since here two consecutive reset events are not allowed, but they must be necessarily separated 
by at least one flight. This has the advantage that the dynamics in the limit r —>■ 1 will be nontrivial since there the 
walker will successively carry out one-flight excursions separated by reset events. Instead, if consecutive resets were 
allowed as in Q, then for r — > 1 the particle would be kept permanently at xg (which looks somewhat unrealistic) 
and the MFPT would diverge in that limit. 

Alternatively, we will compare this with a second mechanism (which has been already used in i) in which the 
statistics of resetting is independent of jumps. In this case flight times are again distributed according to g}(f), but 
the successive times at which a reset occurs are also i.i.d. variables which follow its own pdf 9(f) (independent from 
g^(f)). By introducing this second mechanism (termed as resetting independent of motion) we will be able to check 
if subordination to the motion process has any influence on how resetting affects the values of the MFPT and the 
dynamics of optimal search. 

So that, the evolution of our random walkers will be completely determined for a particular choice of vg and (p(t) (also 
9(t), in the case of motion-independent resetting). To compute the MFPT we will determine first the rate q(t;xg,Xg) 
at which crossings through the origin occur at time t for a free random walk starting from xg. The parameter Xg 
gives the position of the walker immediately after it is reset (for the moment we will set Xg as independent of xg for 
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convenience, while in the end we are mainly interested in the case Xq = xq). Then we can apply the renewal property 

UM 


q{t;xo,XQ) = f{t;xo,Xo) 


Xo, Xo)q{t - t'; 0, Xq) 


( 2 ) 


where f{t]XQ) is the corresponding first-passage distribution through the origin. So trajectories contributing to 
q(t;xo,XQ) are divided into those which cross the origin at t for the first time (first term on the rhs) and those which 
did it for the first time at a previous time t' (last term on the rhs). The connection of the rate function to the free 
propagator, for the case of velocity models with constant speed, is simply given by 


, _ J uoP(0,t;a:o,Xo), it xq ^ 0 

q[T, xo, Xo j I ^ _ 1 if :co = 0 


(3) 


where the additional term —^6{t) has been introduced to neglect crossings at f = 0 by explicitly imposing the condition 
q( 0 ; 0 ,Xo)= 0 . 


The MFPT can be obtained easily now by using for instance Laplace transform techniques [2^ |26|. So, if g{s) 


represents the Laplace transform of an arbitrary time-dependent function g{t)^ then one obtains from ([2]) the expression 
for the MFPT 


(T) 


— lim 

s —>-0 


d/(s;xo,Xo) 


= — lim — 


q(s;xo,Xo) 


— lim 
s->-o as 


ds s-s-o ds [ 1 -I- g(s; xq = 0 , Xq) 

uoP(x = 0 , s;xo,Xp) 


1/2 -I- VoP{x = 0, s; Xo = 0, Xq) 


(4) 


which implicitly provides also the expression for the first-passage distribution /(s; xo, Xq) = q{s; xq, Xq)/ (l + q(s; 0, Xq)) 
in the Laplace space. 

Note that renewal properties like that in ([5]) are strictly valid only for Markovian processes (although they can often 
be used as a reasonably good approximation for more sophisticated cases, specially regarding MFPT computation 
[ 2 ^). So that, in the following we will focus for simplicity on situations which either are Markovian or, equivalently, 
admit a Markovian embedding. 


Resetting subordinated to motion 


The velocity model driven by the process © can be conveniently described through a Continuous-Time Random 
Walk (CTRW) scheme @,111. For this we introduce the probability density ji(x,t;xo,XQ) for the particles starting 
a flight from x at time t. This allows us to write the mesoscopic balance equation 


j(x,f;xo,Xo) = (5(f)(5(x - Xo) 


nt poo 

+ (1 — r) / dt' / dx'dt(x',t')j(x — x',t — t';Xo,X q) 
Jo J — 00 

pt poo 

-f r6{x — Xq) / dt' / dx(p{t')j{x,t — t']Xo,XQ). 

Jo J — 00 


(5) 


The second and third terms on the rhs of account for flights starting after another flight or after a reset, respectively. 
This, together with the initial condition term S(t)6{x — Xq) prevents the particles from being reset until having 
completed at least its first flight. In the expression ([S]) we have introduced the joint probability 'i>{x,t) of performing 
a jump of length x with constant speed vq during time t. Using that flights are done at speed vq or —vq, each with 
the same probability, these pdf’s can be easily related to the pdf of flight times through 


4'(x,f) 



vot) -I- i5(x -I- uot)] ip{t) 



( 6 ) 
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The free propagator for this scheme can be defined as the probability density of particles that are located at point 
X at time t, i.e, 


where we have defined 


nt nOC 

P{x,t;xo,XQ) = / dt' / dx'(l){x',t')j{x — x',t — t'; Xo,Xq)^ 

Jo J — oo 

1 

4>[x,t) = - [(5(a; — not) + (5(a: + Wot)] / dt'ipit'). 

^ Jt 


(7) 


( 8 ) 


This last function gives us the probability that a single flight has not finished yet after having travelled during a 
time t and having covered (either to left or right) a distance not. To find an explicit expression for the free propagator 
we introduce the Fourier-Laplace transform of an arbitrary function g{x,t) as 


g{k, s)= dte 
Jo 


dxe * ^g{x,t). 


(9) 


In the Laplace-Fourier space the expressions ([5][8l) become simplified by virtue of the space and time convolution 
theorems. For example, transforming by Fourier-Laplace Eq. ([5]) we obtain 


j{k,s-,xo,XQ) = e + {1 - r)ii{k,s)j{k,s-,xo,xl) 
+ re~^'"^°j{k = 0,s;xo,xl)if(s) 


( 10 ) 


where the hat and tilde symbols means Fourier-Laplace and Laplace transforms respectively, and k and s the corre¬ 
sponding Fourier and Laplace arguments respectively. Setting k = 0 into Eq. m we can solve for j(k = 0, s; xq, Xq), 
which yields 


1 


j{k = 0,s;xo,xo) = -- 

1 - tp{s) 


( 11 ) 


Finally, by transforming Eq. © by Fourier-Laplace we obtain P{k, s;xo,Xq) = 4){k,s)j(k = 0, s;xo,a;Q), which can 
be combined with Eq. OT to get a closed expression for P{k,s]Xo,XQ) 


P(/c,s;xo,Xq) = 


!j){k,s) 


1 — (1 — r)4'(fc, s) 


g—ifcxo _j_ Xip^s) ^—ikXn 


1 - V>{s) 


where 'I'(fc, s) and (/)(fc, s) are given by 


and 


4f(fc, s) = - [(^(s — ikvo) + (p{s + ikvo)] 


s) = - [(^*(s - ikvo) -I- (^*(s -I- ikvo)] . 


( 12 ) 


(13) 


(14) 


We will focus now for simplicity in the case where the dynamics of flights is Markovian. This corresponds to the 
case when the flight time distribution is exponential, i.e. 


ip{t) = Ae 

Replacing this expression into Eq. m we find after some algebra 

(s -b +r\{l + I) 


P{k,s;xo,xl) = 

and after inverting the Fourier transform 

P{x,s;xo,Xo) = 


(s + A)(s + rX) -b Vgk'^ 


(15) 


(16) 


1 I S + X 


2no \ s + rX 


2sno \ s + rX 


( 17 ) 
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Inserting this expression into (I3B1) we will finally obtain the expression for the first-passage time distribution and 
the MFPT. For simplicity we just show the result for the case of interest Xq = Xq, which read 


/(s : Xo,a;o) = 


• rX 


7-;^ _(_ S g\/(s-|-A)(s+rA)|a:o|/«o 


(18) 


and 


Ar (e-^V^I“olAo 

If we compare this last expression with the results obtained in for the case of Levy flights we observe that the 
scaling hmr_>o(T’) is kept. However, the universal result for xq —t 0 reported there (see Eq. (8) in Reference 

i), i.e. lima;Q_>o(T’)(\/r — r) is not recovered. This is a consequence of the rule ([T]) we have imposed; by forcing 
the walker to make at least one flight before being reset we have introduced some sort of memory in the system (in 
the sense that not all events are equivalent) so in that case the Sparre Andersen theorem does not necessarily hold. 
The different dynamics in both cases is clear by noting that the expression {y/r — r)~^ leads to the divergence of the 
MFPT as r —>■ 1, as we have mentioned above. Instead, it is easy to see from (HU) that the MFPT for our case does 
not diverge in that limit. 


Resetting independent of motion 


If the resetting mechanism follows its own time dynamics (according to the time pdf 0{t), as described above) then 
it is not convenient to use a Continuous-Time Random Walk scheme as done in the previous section to derive the free 
propagator of the walk. Instead, we will use the analogy between random walks with resetting and mortal random 
walks. Following the nomenclature of several recent works mortal random walks are defined as walks subject 

to a mortality function such that after a random time, governed by a given pdf, the particle will disappear (it is, it 
will ’die’). So, if we define in our random walks with resetting an excursion as the action that goes from t = 0 to 
the first reset event or, alternatively, from the i-th to the i + 1-th reset event, then every one of these excursions can 
be interpreted as an independent mortal random walk. Consequently, the random walk with resetting consists of the 
time convolution of successive mortal random walks governed by the ’mortality’ function 9{t). 

We will focus again for simplicity in the Markovian case, so we choose for the resetting mechanism 9{t) = 
with uJm being the frequency of resets (or, alternatively, the ’mortality’ rate). In addition we consider ip{t) = Ae“^*, 
as done in the previous Section. For this specific choice, the free propagator of a mortal random walk in the Laplace 
space for an isotropic velocity model has been derived before and reads 


Pm{x,S]XQ) 


2vo 


^ 1 S+UJm+'^ 


■ LOr) 


( 20 ) 


where the subindex m stands for mortal. 

Using the ideas discussed above and denoting the time convolution operator as *, we can write the free propagator 
for the random-walk with resetting as 


P{x,t-,xo,xl) 


= Pm{x,t-,Xo) + Pm{x,t;xl) * ip{t) 
+ Pmix, t; xl) * (fit) * ip(t) + ... 


or, transforming to the Laplace space, 


P{x,s] xo,Xq) = Pm{x,s;xo) + 


Pm{x,S-,Xo)ip{s) 
1 - V3(s) 


( 21 ) 


( 22 ) 


The first term on the Ihs of (ED) accounts for the case where no resetting has occurred yet since t = 0, the second 
term stands for the case where only one resetting to the position Xq has occurred up to now, and so on. 

Now, following the same procedure as above, we can insert Equations (1^ and (1^ into (|3]) to provide an expression 
for the first-passage distribution 


f{s;xo,XQ) 


■ 


+ 5(1 + 




(23) 
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and for its corresponding MFPT 


(T) 


1 


1 + 


OJm 





(24) 


One remarkable point is that one can recover from this result in the diffusive limit (A —>■ oo, uq —>■ oo with Uq/A —^ D) 
the expression 

(T) = — (eV^^rn/DM _ (25) 

tOm ^ ^ 

which was derived for a Brownian walker with resetting in [l[ . Likewise, it is easy to check from (IMl) that lim^j^ ^o{T) = 
oo and lim,^^_>oo(’r) = oo, which guarantees the existence of an optimum resetting rate, except for the case Xq = 0, 
where the MFPT decays monotonically with and so the minimum MFPT is obtained for oo. This 

qualitatively coincides with the result found for the case of resetting subordinated to jumps, which shows that in the 
limit xo —>■ 0 the optimal strategy is always resetting as much or as faster as possible. A more complete discussion on 
this point is provided in the next section. 

_I /o 

Again, in the limit oJm —>■ 0 we obtain a scaling (T) ^ Um ^ , which really seems to represent a universal behaviour 
independent of the mechanism used for resetting. This universal feature is actually reminiscent of the scaling ~ 
for the revisit times to a given position of random walkers which have jump pdf’s with finite moments. The first-passage 
statistics is expected to be asymptotically governed by the probability to be brought back (through a reset) to the 
initial position after very long excursions. Specifically, MFPT should be inversely proportional to that probability, 
which leads to (T) ^ with (t) the mean time for resetting. Consequently, the scaling (T) ^ follows 

straightforward. 


OPTIMAL SEARCH TIMES IN INFINITE DOMAINS 

Next we will use the results m and (1241) just derived to analyze the existence of an optimal point in the phase 
space (r. A) (or (wmjA)) such that it minimizes (T) as a function of xq. Since we have deliberately chosen a very 
simple scenario in which both processes (resetting and motion) are characterized by one single scale (together with 
the spatial scale of the problem, xq), one should expect that the optimal strategy is relatively trivial, in comparison 
with more complex multi-scale or scale-free situations as in Levy flights, which surprisingly lead to a first-order 
transition in the phase space for the optimum (T) Q. However, it is still of interest to check whether the different 
resetting mechanisms proposed will lead to a very different behavior. As we will show in the following Section, the 
dynamics becomes specially interesting when finite domains are considered (so a new spatial scale is introduced) with 
the emergence of phase transitions in the phase space (r. A) or (ujm, A). 

We start by analyzing the result (I19p . By differentiating with respect to A one finds that a minimum is found for 

A - 

" V^\xo\ 

where W() denotes the Lambert-W function. One can check that Aopt is always positive for any value in the interval 
0 < r < 1. Replacing its expression into m the global minimum is always found for r = 1, which corresponds to 
^opt = yco/la^ol: where we define y = (1 -|- W(—0.5/e)) ~ 0.768. So that, the global strategy in this case corresponds 
trivially to using resetting as much as possible while adapting the flights durations to the spatial scale xq. 

Next we go for the case of resetting independent of motion. Since one can easily check from (1241) that the derivative 
of the MFPT with respect to A is always positive, this immediately leads to the conclusion that the global optimum 
must be at A = 0. Differentiation with respect to LOm then allows one to hnd that the optimum satisfies (<j-'m)opt = 
yxo/lxol ~ 0.768uo/|a^o|. It is found then that the optimal resetting rate uJm coincides with the case of resetting 
subordinated to jumps. This actually shows that the optimal strategy for the random walker in both cases is exactly 
the same and consists of single-flight excursions interrupted by resetting at the rate yuo/|2;o|- These results are 
actually very intuitive: resetting will be in general a more convenient mechanism for going back to the initial position 
than turning the direction of motion (at least while resetting is considered to be an instantaneous process), and 
consequently optimizing the MFPT should rely just on appropriately adjusting resetting to each particular situation. 

We stress that additional mechanisms for resetting have been tested in which the same optimal behaviour has been 
obtained, again with the same optimal resetting rate yco/lxol. So, we can conclude from the above that this is a 
universal feature for all random walk models with constant velocity vq provided that resetting is Markovian. 


1 -k W 


-1 


e(l -I- ^/r) J 


(26) 
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OPTIMAL SEARCH TIMES IN FINITE DOMAINS 

We consider now that the random walkers move in the interval (0, L) with periodic boundary conditions. Note that 
the procedure to obtain the MFPT described in ([3]l4]) is still completely valid, provided that now the free propagator 
P{x,t;xQ,XQ) is replaced there by the propagator Pi(a;, t; ccq, Xq) for a finite domain of size L. Transforming from one 
to another is relatively simple just by explicitly imposing the periodic condition at the boundaries [i^ . 

OO 

PLix,t]Xo,XQ) = ^ P{x + mL,t]Xo,xl), (27) 

m— — oo 


which is valid in the region of interest x G (0, L). 

Introducing in (l27l) the expressions for the free propagators we have already derived in previous Sections one can 
obtain again the corresponding MFPT. Since these calculations do not provide any additional insight into the problem 
we will just reproduce here the final expressions obtained. For the cases of subordinated and independent resetting 
with Xq = xo the MFPT reads, respectively, 



1 + + (1 - 

Q — \^xo/vo _|_ g—Av^(L—xo)/«o 


- 1 


(28) 


and 


(T) 


(l + (l - 

3--y/aJm(‘^m+A)xo/'i'0 _j_ g-\/l^m(<^m+A)(L-Xo)/t)o 


- 1 


(29) 


from which (I19p and (I24|) can be recovered, respectively, in the limit L —>■ oo. 

The main focus of interest now is in understanding how the global optima found in the previous Section for infinite 
domains changes as the domain size gets reduced. For this, note that in the case where the spatial scales xq and L 
are of the same order then both resetting and changes of direction will be detrimental, since coming back to Xq is 
only a good strategy provided there is a target close enough to there, and this will only be the case for xq L (or, 
equivalently, |L — xo| L). Hence, as long as L decreases and that condition is not fulhlled, we should find that the 
optimal point in the phase space (r, A), or (wm, A), will move towards (0, 0). 

Interestingly, what we obtain from both (l28l) and (|2^ is that the transition from one regime to the other is not 
smooth but follows a first-order phase transition. Figure [T] provides plots of (T) as a function of r and A for resetting 
subordinated to motion. There it can be seen how the point (0,0) represents a global minimum for xq/L large 
(panel (d)), but when this parameter is small an intermediate minimum appears and eventually it becomes the global 
minimum (see panels (a) and (b)). This picture can be completed with the information provided in Figures [5] and 
|3] (which correspond to the case subordinated to motion and independent of motion, respectively). There we have 
computed the optimal value (of A or Wm, respectively) as a function of the order parameter xq/L, and we have marked 
with a vertical dotted line the point at which the phase transition occurs. The scaling Aopt ~ xff^ and (a;m)opt ^ 
can be also observed there. 

The critical point can be easily approximated as follows. First, one should see that the MFPT at the point (0,0) 
reads (T) = L/2vq, as can be computed from (l28l) or (|2^ : this coincides with the MFPT one would expect for a 
ballistic particle moving in the periodic domain (0, L). Then we should determine when the intermediate minimum 
observed in Figured] takes a value equal to LJ^vq. For doing this, we note that the coordinates of that point can be 
well approximated by the optimum point we have found in the previous Section (for the case of infinite domains); the 
justification for this is that we still expect the condition xq <C L to be valid close to the critical point. Then we will 
use the points (r, A) = (I,xz;o/xo) for the subordinated case, and (wmjA) = (xxo/xo,0) for the independent case. If 
we replace this into (IT^ and (l24l) we obtain 

(t)„= 2(2^31) (30) 

as the critical point at which the transition occurs for both mechanisms. This result coincides well with that reported 
directly from Eqs. (I28I29|) in Figures [2] and [3|(see insets there). Again, we stress that this result will be valid for any 
Markovian resetting mechanism provided that random walks are isotropic and the walkers move at a constant speed. 
















Optimal times with delayed resetting 


We find then that our walkers with constant speed must adopt a very extreme strategy in order to optimize its 
search times, as they must always interrupt excursions with a reset just after the first flight is completed. One may 
think that this is simply because resetting is considered here to be an instantaneous process, so it is costless in terms 
of time. So, we could introduce a delay such that when a reset is carried out the particle must wait an average time r 
before starting a new excursion. This would punish resets so a different (maybe richer) dynamics could be expected. 

We have explored this possibility both in the resetting subordinated to motion and independent to motion, but just 
show the results of the former for simplicity. In particular, we can replace the last term in Equation ([5]) by 

pt nOO 

rS{x — Xq) / dt' / dxip{t')j{x,t — t' — t;xq,Xq), (31) 

Jo J — oo 

so we explicitly introduce the delay t in the resetting term. With this change now the free propagator in the infinite 
domain (Eq. (IT7l) l reads 


P{x,s;xo,xl) 


-■^/(5+A)(s+rA)|x-3;o|/^;o 

2 V S + r A 

_ ^ _ I s + X 

2 (se"^® + rA(e'^® — 1)) V s + rA^ 


s+A)(s+rA)|x—Xpl/iiQ 


(32) 


Using again Eq. (EH) and repeating the same procedure as above for obtaining the MEET, the result reached for 
Xq = Xq is 


(T) = (l + Arr)(ro), (33) 

where (Tg) is the MEET in the absence of the delay (given by the expression (12^ 1. Actually, it is possible to extend 
this result (not shown here) to check that (1551) is also valid when the time delay is a random variable that follows any 
pdf with average t. Likewise, the case of resetting independent of motion leads similarly to 

(T) = (1 + Ar)(ro), (34) 

Since the delay only introduces a multiplicative factor in the computation of the MEET, the effect of this on 
the optimal strategy is actually minor. One can check (Eigure 01 triangles) that the optimal point in the phase 
space (r. A) for resetting subordinated to motion decreases when the delay r increases, and eventually for a critical 
value of T the phase transition disappears and the trivial point (0,0) remains as the unique global optima. For the 
resetting independent of motion (Figure 01 circles) the delay does not introduce any change in the dynamics of the 
optimal MEET, since the parameter uJm does not appear in the prefactor in (155)) . This reflects that the equivalence 
between both reset mechanisms studied here break down when additional (non-Markovian) effects as a delay effect 
are introduced. 


DISCUSSION: LINKS BETWEEN RESETTING AND MULTI-SCALE WALKS 

First-passage time distributions of simple random-walks in infinite media typically decay very slowly following a 
power-law function or similar, the most famous case being the scaling ^ found for isotropic Markovian walks, a 
result known as the Sparre Andersen theorem. Instead, resetting destroys this scaling by introducing memory effects 
in the process and leads eventually to an exponential decay of the first-passage pdf, as typically found in finite media 
(albeit one must be aware that both situations are of a different nature; while the latter reaches an equilibrium state 
the former asymptotically approaches a nonequilibrium stationary state). As a consequence of this, we have showed 
here that an optimal resetting rate (in terms of minimizing the MEET through the origin) does appear. From that 
point of view, our work extends the case of diffussive movement with resetting which was already explored in references 
[i|, 0 to the case of walkers moving at constant speed vg and carrying out exponentially distributed flights. It is so 
illustrative to recall the optimum resetting rate found in [l|, which in our notation corresponds to 2.538uQ/Aa:g. Here 
we have checked that the case of exponential flights leads in all cases to an optimal rate 0.768uo/|a;o|. The different 
scaling found in the two cases comes from the fact that our random walkers follow a velocity model, so characteristic 









9 


times for the motion process over a distance x are typically x/vq (in comparison to Brownian motion where typical 
times read x^ jD = Xx^/vq). 

The case of finite domains we have also explored (now with xg defined as positive) shows that when we introduce 
the additional spatial scale L, then Xq/L plays the role of an order parameter such that a first-order phase transition 
for the optimal resetting rate is found. So, for xq/L small we still can approximate the optimum rate by 0.768uo/xo, 
while above the critical point xq/L ~ 0.1159 the optimum resetting rate is trivially zero. It is interesting to note 
that an analogous behavior has been found for the case of two-scale random walks (without resetting) in finite media 
[ 2 ^ . If the walker is allowed to move according to two characteristic movement scales (say Ai and A 2 ), then we 
also find a critical value of xq/L ~ 0.105 such that above this value the global optima corresponds to Ai = A 2 = 0. 
Instead, for the region under the critical value the best strategy corresponds to Ai = 0 and A 2 ~ O.Sxo/xq. So that, 
we observe that the role played by this second movement scale A 2 is similar to the effect that resetting has (although 
both mechanisms have fundamental differences, since resetting clearly introduces a bias in motion towards the initial 
point while the two-scale walk is completely isotropic). These results reinforce the idea discussed in previous works 
on search theory 3, 24, ^ that optimizing search efficiency of random walks (in terms of minimizing their MFPT) 
necessarily implies an appropriate combination of different (at least two) motion scales, which in the limit case leads 
to the use of free-scale, i.e. Levy, strategies. Smaller scales will then be used to efficiently exploit closer regions 
when this contains targets, while larger scales will be used to explore further regions. While a more generalized and 
fundamental study of the link between multiscale motion and MFPT optimization is still elusive, the present work 
provides so novel and significant evidence in this line. 
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FIG. 1: MFPT for resetting subordinated to motion as a function of the resetting probability r and the average flight length A 
for different values of the spatial scale: Xo/L = 0.005 (a), 0.025 (b), 0.1 (c), 0.25 (d). The legends on the right of each panel 
provide information about the order of magnitude of typical values of the MFPT. Arbitrary values uo = 1 and L = 200 have 
been used in all the cases. 



xo/L 

FIG. 2: Optimal value of the rate A for the case of resetting subordinated to motion. The vertical dotted line represents the 
critical value for which the global optima becomes (r, A) = (0,0). The behaviour close to the critical value is shown in the inset. 
Arbitrary values uo = 1 and L = 200 have been used in all the cases. 
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FIG. 3: Optimal value of the rate A for the case of resetting independent of motion. The vertical dotted line represents the 
critical value for which the global optima becomes (aJ„i,A) = (0,0). The behaviour close to the critical value is shown in the 
inset. Arbitrary values wo = 1 and L — 200 have been used in all the cases. 



T 


FIG. 4: Critical value oi xq/L as a. function of the delay r introduced after the reset. Results are shown both for the case 
of resetting subordinated to motion (triangles) and resetting independent of motion (circles) to observe the different dynamics 
that the delay induces in each case. Arbitrary values va = 1 and L = 200 have been used in all the cases. 






